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Predicting the incidence of planets and debris discs as a 
function of stellar mass 
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ABSTRACT 

The mass of solids in a young circumstellar disc may be the key factor in its 
efficiency in building planetesimals and planetary cores, and dust observed around 
young T Tauri and Herbig Ae stars can be used as a proxy for this initial solid 
content. The dust-mass distributions are taken from recent millimetre- wavelength data 
and fitted using survival analysis to take into account upper limits, and threshold 
disc-masses for building planets and belts of comets are estimated. Amongst A-stars, 
20 % gas giant and 55 % debris disc systems are predicted, in good agreement with 
observations. For M-stars, the predicted and observed planet-frequencies agree at ~ 
2 — 3 %, and this low incidence is explained by a lack of massive discs. However, debris 
is predicted around ~ 14 % of M-stars, while only ~ 2 % such systems have so far 
been found. This suggests that deeper searches such as with Herschel and SCUBA-2 
may find a cold disc population previously missed around these low-luminosity stars. 
Also, an estimate of the efficiency of building millimetre-detected dust into planetary 
cores suggests that about a third of M-stars could host an Earth-mass planet - but as 
the dust is spread over large disc areas, such planets may orbit far from the star. 
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1 INTRODUCTION 

The incidence of extrasolar planets is known to be a function 
of both stellar mass and metallicity. If the present-day metal- 
licity of the star reflects that of its circumstellar material at 
early times, then metal-rich discs would be expected to build 
planetary cores more readily by grain coagulation. The stel- 
lar mass may enter if more massive stars have higher-mass 
discs, as well as in factors such as shorter dynamical times 
for grains to collide. In the simplest approach, the mass of 
solids in a disc, i.e. the total mass M multiplied by the solid 
fraction Z, may be the relevant threshold for a successful 
outcome in forming plan ets or planetesimals (|Wvatt et alj 
120071 : iGreaves et alj|2007l ). These two quantities M,Z can be 
readily estimated using proxies. Millim etre- wavelength emis- 
sion fr om dust is rather optically thin l|Andrews fc Williams! 
l2007al ) and the derived disc masses in present-day star for- 
mation regions are canonically multiplied by a factor « 100 
to include the gas content inherited from the interstellar 
medium. These masses are proxies in the sense that mate- 
rial already converted into larger bodies (from boulders up 
to planets) has little emitting area, and so potentially only 
some of the orbiting material is represented. A good proxy 
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for the refractory proportion of material in discs of non-solar 
metallicity is the logarithmic [Fe/H] of the host star. 

IGreaves et al.l l|2007l ) showed that such a simple empir- 
ical model can reproduce the properties of both planets and 
debris discs around Sun-like stars. The presence of debris 
indicates collisions between planetesimals, such as icy outer- 
system comets, treated in the model as a less successful mode 
of planet formation ending only in small bodies. Both the 
frequencies and metallicity-dependencies of these two out- 
comes were well matched. In particular, forming planets re- 
quires a high solid mass to build the core, and so high metal- 
fractions are advantageou s, explaining the s trong metallic- 
ity trend in the host stars (|Wvatt et al.ll2007h . On the other 
hand, forming comet belts requires only modest amounts 
of material, and many discs meet this thr eshold regard- 
less o f metallicity, so it enters only weakly (IGreaves et alj 
120071). More complex approaches such as the simulations of 
IJohansen et al.l (|2009l ) have confirmed that the outcomes of 
pla netesimal formation sh ould depend strongly on metallic- 
ity. iRavmond et alj l|2007l ) have used a model of disc-masses 
in solids to investigate the formation of terrestrial planets, 
finding in particular that habitable examples could be less 
likely for low-mass stars if they had insubstantial discs. A re- 
cent summary of our understanding of p lanetesimal growth 
is presented bv lChiang fc Youdinl |2010) 
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Figure 1. Distribution of disc masses (in gas plus dust) for A, 
FGK and M stars. Disc detections and upper limits are shown 
by filled and unfilled bars respectively; the red lines shows the 
fitted population from survival analysis (see text). X-axis labels 
are for the bin-centres. For reference, a solar mass is just over 1000 
Jupiter masses; the lowest bin comprises discs of < 1.3 Mj up ;t er 
or <4 M-Earth in solids at solar metallicity. 



Here our empirical model is expanded to test the planet 
and debris frequency expected for stars of higher and lower 
mass, namely A and M dwarfs. Their planetary systems 
are now being discovered by radial velocity surveys (of sub- 
giant descendents in the case of A-stars), and their debris 
disc s seen as far-in f rared excesses e.g. b y the Spitzer satel- 
lite l|Su et al.ll2006l ; lGautier et alj|2007l ). Here the M,Z dis- 
tributions are assessed, and while the [Fe/H] ranges among 
the stars are all similar, the disc-mass distributions differ 
markedly, leading to strong changes in incidence with stel- 
lar mass. 



2 DATA 

Metallicity data were first compared for different stellar 
types. Subgiants that are the descendents of A-stars of ~ 1.2- 
2.5 solar masses can be identified by absolute magn itude 
My of 0.5-3.5 and B-V of 0.55-1.0 d Johnson et al.ll2006l). o r 
similarly effective temperature 500 0-6300 K jFlowerl 1996T) . 
The SPOCS study of nearby stars (|Valenti fc Fischerll2005T l 
yields approximately 100 such subgiants with [Fe/H] mea- 
surements. A (log-)normal fit to these values indicates mean 
and standard deviation /i, a of +0.02, 0.23. This is very 
similar to values for nea rby Sun-like dwar f s, wit h fi, a of 
-0.01,0.29 in the data of IValenti fc Fischerl (120051 ). For M- 
stars, the derivation of metallic ities based on photom etry 
has recently be re-computed by I Johnson et alj (|2009h . Us- 
ing data they present for M-stars within 20 pc, and their 
relation for [Fe/H] derived from K-magnitude height above 
the main-sequence, fi, a are found to be -0.06,0.24. These 
distribution^ are so similar that they can not significantly 
affect planet-building potential around different classes of 
star. 

The disc mass distributions were then constructed and 
found to be quite different for A versus FGK and M stars 
(Figur e 1). The masses for the A-stars are from lNatta et all 
(2004), in a millimetre-based compilation for young Herbig 
Ae objectfl Their plot of M-disc/M-star versus M s tar was 
used to derive Mdisc for M star of 1.2-2.5 solar masses. Since 
11/45 of the disc-mass measurements were upper limits, a 
survival analysis was used to find the unde rlying distribu- 
tion, implemented with the asurv package l|Lavallev et al.l 
1092 .. The resulting disc-mass distribution appears rather 
log-normal in form, and is centred at a mean of 25 Jupiter 
masses of gas and dust. For comparison, this is just above 
the Minimum Mass Solar Nebula of ~20 Jupiter masses 
(Davis 2005) needed to supply the cores of all the Sun's 
planets. For FGK and M starjfl, the millim e tre-su r veys o f 
Tau/Aur and Oph by lAndrews fc Williams! l|2005l . l2007bh 



1 ISchlaufman fc Lauehim] feOlOT) argue that the mean [Fe/H] for 
FGK and M stars are about 0.1 dex lower than found here; such 
small shifts have negligib le effect on the M. n udx d istribution. 

2 The data presented by [Alonso-Albi et al.l J2009I) give a similar 
mass-distribution but with fewer HAe stars. 

3 Pre-main sequence M-stars have rather vertical Hayashi tracks 
with modest change in effective temperature, so young and old 
objects of M-type will have similar stellar masses. Other pre-main 
sequence stars evolve towards hotter spectral types, but since 
most of the observed T Tauri objects have K-classifications, they 
will still be within roughly 'Sun-like' classes at later times. 
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Tabl e 1. Predicted and observed percentages of giant planets and debris discs. The planet data are f r om the uniform-detectability sample 
of Ijohn son et all ll2010h and debris frequencies are from Spitzer 70 urn surveys l|Gautier e t al. 2007; Trilling et al,| [2008: Bcichma n et al.1 
2006; Su et all2006h . limited to data with signal-to- noise > 3. The debris frequency observed for FGK stars is used to make t he predictions 
for o ther spectral types. The upper limit for M-stars is from <1/15 objects detected by Spitzer at 70 umin ^ 3<r data l lGautier et al.l 
2007), while the bracketed value is from 1/50 detections in the millimetre survey of iLestrade et al.l J2009T) . Poisson errors are based 
on the number of planet and debris systems detected, and the numbers of discs falling within the input mass-bounds (see text). The 
parameters of Equation 1 are list ed in the four colum ns after the spectral types. The A-star mass marked with a *-symbol is actually the 
se cond most-massive value from lNatta et al as another disc with large distance and uncertain mass was removed. The sample 

of lAlonso-AIbi et al. (2009) gives a comparable highest mass of 225 Jupiter masses among HAe stars within about 500 pc. 
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were used, c ompr ising similar-aged stars to the HAe's of 
iNatta et al.l l|2004 ). Since these authors used a conversion 
of flux t o dust mass yieldin g values a factor of three lower 
than in INatta et all (|2004h . the masses have been shifted 
upwards by one bin (51og-M = 0.5) in Figure 1. The dif- 
ference arises from using an opacity at 1 mm wavelength of 
0.03 cm /g rat her than the . 01 cm /g found in the detailed 
grain model of INatta et al. (120041) . which also agrees with 
recent literature values (|Drainell2006j ). Survival analysis was 
again used to construct the base disc-mass populations for 
FGK and M stars, but more of the points are now upper 
limits, as the discs are less massive than for the HAe stars . 
There were 57 limits out of 111 measurements for the FGK 
dwarfs, and 57/106 for the M-stars. The means are 4 and 
3 Jupiter masses respectively. There is also a large tail of 
insubstantial discs in both cases, which combined with the 
lower means, suggests that these discs will have much less 
planet-forming capability than those of A-stars. The low- 
mass tails include many Class III T Tauri sta rs, which may 
lose their discs early on jLuhman et alj|20ld ). All this disc- 
mass distributions are very broad, whereas the metallicity 
distributions would cover only about one bin on the scales 
of Figure 1. 



3 RESULTS AND PREDICTIONS 
3.1 Planets 

Since the metallicity distributions are all narrow and simi- 
lar, the differences in Mdisc should be the dominant factor, 
and the r elevant thresho l ds for planet formation are now 
estimated. iGreaves et al.l |2007l ) argued that if a minimum 
solid mass is needed to grow a planetary core, then there will 
be observed examples of giant planets arising from a range 
of discs with this same solid-mass, but ranging between 
massive/low- metal and lower-mass/high-metal systems. The 
lowest-metal system observed with a planet arises when the 
solids criterion was met only because the disc was at the 
extreme high-end of total mass. The solid mass can then be 
expressed using the proxy data as 



Tirrmn 

solids 



in /100. 



(1) 



The resulting parameters for all the spectral types are 
listed in Table 1. The input data on the disc masses are 



from the references in the previous section, and the metal- 
licities of planet-hosts are from http://exoplanet.eu/ The 
subsequent columns list the disc-mass thresholds needed to 
form planets, and the predicted and observed frequencies. 
For FGK stars, the predicted incidence of around 8 % is 
in excellent agreeme nt with the 8.5 % foun d in the system- 
atic planet sample of Ijohnson et aL I (|2010l). For A-stars the 
planet agreement is also very good, with around 20 % both 
predicted and observed. For M-stars, the incidence is 2-3 %, 
and the model correctly predicts a low value because few of 
the discs are massive (Figure 1). 

The solid-mass thresholds needed for forming giant 
planets are here found to be 200, 400 and 650 for A, 
FGK and M spectral types respectively (Table 1). This trend 
qualitatively matches the expectation that grains will coag- 
ulate more readily in discs around more massive stars, as 
the dynamical (orbital) times are shorter. 

Subsets of different types of planetary systems are not 
considered in detail here. Rocky 'super-Earth' planets form 
a small part of the FGK-hosted population (e.g. < 10 M® 
bodies contribute around 0.5 % incidence), and the require- 
ments for building such bodies are considered briefly below. 
Also, there are likely to be system s with more gas g iants 
orbiting further out - for example, iGould et al.l (|2010l ) es- 
timate ~ 1/6 for the frequency of two-giant systems, from 
microlensing results - but these are not included here. The 
analysis is however self-consistent in using the lowest-metal 
radial-velocity detection to predict the planet frequency in 
this type of survey. 



3.2 Debris discs 

For the less-successful outcome of building planetesi- 
mals, the thre s hold can not be straightforwardly derived. 
IGreaves et al.l (|2007l ) used the metallicity range of FGK 
debris-hosts to estimate the mass-threshold for building 
comet belts, but there is little similar data for A and M- 
stars. Hence the simplest procedure is to infer the thresh- 
old from the data. In deep Spitzer surveys, approximately 
19 % of FGK stars hav e debris discs detected at 70 [im 
l|Greaves fc Wv att 2010), associated with belts of colliding 
planetesimals. Associating this outcome with discs of lower- 
mass than those that build planets, this next 19 % of the 
Mdisc distribution corresponds to >50 Earth masses in solids 
or >15 Jupiter masses in total at solar metallicity. This is 
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eight times lower than the value needed for building plan- 
ets. Making the minimal assumption that the efficiency of 
the planetesimal-building process is similar in all systems, a 
similarly reduced value can be applied for the other stellar 
masses. Discs of at least 25 and 80 Earth masses in solids 
would then be needed for A and M stars to have debris, 
respectively. 

The model results (Table f) show a predicted debris 
incidence for A-stars that is slightly low, but only at the 
-lcr bound (error on observed— predicted of 15 % in quadra- 
ture). For M-stars, there is however an interesting discrep- 
ancy between the debris prediction of around 14 % and the 
observed incidence of as little as 2 % (Table 1). A few other 
M-host debris discs have been discovered, but all around 
stars much younger than typical main-sequence lifetime s 
|Liu et alJl2004lForbrich et alJliooll ; iPlavchan et al.ll2009l) . 
It is possible that these M-discs evolve down to lower dusti- 
ness than seen among A or FGK stars (Figure 2), but given 
that FGK and M discs are similar at the T Tauri phase, 
an alternative explanation is favoured. As M-dwarfs are of 
low luminosity, dust particles in thermal equilibrium with 
the star will be cool, compared to a grain on the same orbit 
around AFGK stars. In this case, the discs will have faint 
emission shi fted to long wavelengt hs, and this introduces a 
survey bias l|Matthews et al.ll2007l , e.g.). There may be neg- 
ligible excess above the M-star photospheres at 70 urn, and 
exploratory millimetre surveys may not yet have reached 
deep enough to detect the discs. 

There are some systems known with both gas giants 
and debris discs, but th e statistics suggest t he two phenom- 
ena are not correlated (jBrvden et al.l [2009h : here in fact it 
is assumed that they are separate outcomes. In the uncor- 
rected case for FGK stars, a proportion of 0.08 (planet- 
frequency) x 0.19 (debris frequency) or 1.5 % of stars would 
have planets and dust, so this should be a negligible popu- 
lation. A small number of planets around A-hosts of debris 
discs has also been found in direct imaging surveys. These 
cases are not considered as a separate class here. The sta- 
bility of debris discs over several Gyr, in terms of size and 
m ass and largest bod i es pre sent, has been studied in detail 
bv lHeng fc Tremainel (|2010T ). 

3.3 Planet-building efficiency 

The planetary systems of M-dwarfs are of particular inter- 
est because these stars are so abundant in the Galaxy. The 
large solid-threshold inferred for M-stars to form planets 
(Table 1) suggests the core accretion process may be in- 
efficient. A rough efficiency estimate can be made, assuming 
th e solid content i n the giant planets in the M-star sample 
of IJohnson et all (|2010T ) is of order 20 Msarth- For exam- 
ple, these planets a re of > Msatum and th is is about Sat- 
urn's solid content (|Saumon fc Guillot|[200l l. A lower-mass 
'super-Earth' population is also being discovered of which 
the most massive observed in tran sit has a 20 ME ar th core 
plus RilO-15 % atmospheric mass l|Bean et al.ll200Sl ). If 20 
Earth-masses is a rough boundary above which planets got 
on to become gas giants, then we could equate this outcome 
at ~ 3.3 % frequency with the same fraction of top-end discs, 
which have > 300 M® of dust. In this case, the efficiency of 
incorporating disc solids (or at least the millimetre-detected 
component) is ~ 20/300 or 7 %. 
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Figure 2. Distributions of percentage of stars versus fractional 
dust luminos ities for A and FGK systems with debris detected 
at 24 ,70 um jSu et alj|200ri iTriUing et aljliooj ; iBeichman et al.l 
l200fjt ). The t y pical upper limit for the M-stars observed by 



iGautier et alj I 20071 ) is shown at the top, and applies for dust 
at temperatures of 40-100 K. 



One consequence is that to build an Earth-mass planet, 
15 Mq of dust in the disc would be needed. This occurs 
for around 1/3 of M-stars (Figure 1), so they are not all 
expected to form Earth-analogues. Given that about 3/4 
of main-sequence stars are of M-type, at least a quarter of 
nearby stars could host an Earth-mass planet under this sim- 
ple scaling app roach. However, the disc mat erial spreads far 
from the star (I Andrews fc Williams! l2007bh . so these plan- 
ets may only be able to form in large orbits, particularly as 
the Hill radius for accretion increases with radial distance. 
Microlensing is particularly well suited to discovery of M- 
dwarf systems ho sting such low-mas s planets in orbits out 
to around 10 AU (jGould et al.ll2010h . 

It is beyond the scope of this study to assess how many 
planets might migrate inwards, to the the habitable zone 
around 0.1 AU wher e temperatures are su itable for liquid 
water on the surface. iRavmond et~all (|2007m have simulated 
the formation of terrestrial planets, and find that habitable 
planets of > 0.3 M® (the boundary for plate tectonics and 
atmospheric retention) are much less likely for M-dwarfs 
than G/K stars. Figure 1 also suggests that discs substantial 
enough to form an Earth-mass planet are much more com- 
mon around A stars, although these discs may also form 
giants, with implications fo r the supply of mater ial and sta- 
bility for terrestrial bodies (|Fogg fc Ne lson 2009, e.g.). 



4 CONCLUSIONS 

The disc-mass distributions for young A and M stars have 
been constructed and compared to those for Sun-like stars. 
Inferring the mass-thresholds for forming extrasolar planets 
and debris discs, the higher incidence of both phenomena 
among A-stars or their descendents is well reproduced, as 
is the low incidence of giant planets around M-stars, which 
rarely have substantial discs. However, only about 2 % of 
nearby M-stars have debris discovered so far, and several 
times more are expected. There may be a population of cool 
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debris discs around red dwarfs largely missed in su rveys so 
far. The new SCUBA-2 camera l|Holland et al.ll2006h observ- 
ing at 450 and 850 um will be ideally matche d to cold dust, 
as is the highly sensitive Herschel satellite (|Pilbratt et al.l 
2010), now observing out to 500 urn. A further prediction of 
the empirical calculations is that about a third of M-stars 
have the capability to form a 1 Earth-mass planet - how- 
ever, such planets are likely to form far from the star, as the 
disc material is very spread-out. 
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